setwd("D:/Rev2Plot/")
rt<-read.table("Derived.dat", header=TRUE)
tmc<-read.table("Derived_TMC.dat", header=TRUE)
cmc<-read.table("Derived_CMC.dat", header=TRUE)
pmc<-read.table("Derived_PMC.dat", header=TRUE)
par(mfcol=c(3,1))
oldpar <- par(mar=c(5,5,1,0.8))
hist(rt$Tau,breaks=seq(0,2,by=0.2),main="ALL",xlab=""                                 ,ylab="N")
hist(tmc$Tau,breaks=seq(0,2,by=0.2),main="TMC",xlab=""                                 ,ylab="N")
hist(cmc$Tau,breaks=seq(0,2,by=0.2),main="CMC",xlab=expression(Optical~Depth(phantom()^{13}*C*O)),ylab="N")


